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The surface impedance Z{f) of conventional isotropic materials has been carefully measured for 
frequencies / ranging from IkHz to 3 MHz, allowing a detailed investigation of the depinning tran- 
sition. Our results exhibit the irrelevance of classical ideas on the dynamics of vortex pinning. We 
propose a new picture, where the linear ac response is entirely governed by disordered boundary 
conditions of a rough surface, whereas in the bulk vortices respond freely. The universal law for Z{f) 
thus predicted is in remarkable agreement with experiment, and tentatively applies to microwave 
data in YBaCuO films. 



PACS numbers; 74.60. Ge, 74.25.Nf. 

A perfect sample of a type-II superconductor in the 
vortex (or mixed) state would be transparent to an elec- 
tromagnetic wave at very low frequencies. But defects 
are always present and strongly alter the quasistatic 
and low-frequency response; low frequencies here means 
n = 271 f <C r^d, a so-called "depinning frequency" |^] 
depending on the material and vortex density. It is im- 
portant for applications to know what kind of defects 
can pin vortices, how they hinder small vortex oscilla- 
tions and thereby restrain the penetration of an ac rip- 
ple. In this respect, a study at low levels of excitation 
of both the resistive and inductive part of the surface 
impedance Z{Cl) = R — iX as a function of the frequency 
provides much information about the dynamics of pin- 
ning. It is generally accepted that bulk pinning centers 
limit the quasistatic skin effect to a pinning (or Camp- 
bell's) length Ac ~ 1 — 100 fim, while dissipation is van- 
ishingly small, as observed 0,1). No model however has 
been able to account for variations of Z at higher fre- 
quencies. In particular, as the first increasing of R{f) is 
stronger than expected, the understanding of dissipation 
remains a puzzling problem, including in high Tc mate- 
rials 1). 

Experiments are performed on a series of slabs of cold- 
rolled polycrystalline Pbin and pure single-crystralline 
Nb. The slabs (xy) are immersed in a normal magnetic 
field B; unless specified their thickness 2d is much larger 
than the flux-flow penetration depth Sf (see below). At 
equilibrium, up to the upper critical field Bc2, a regu- 
lar lattice of vortex lines parallel to z is formed, with 
the density n = B/(po, where ipo is the flux quantum. 
Both faces of the slab, z = ±d, are then subjected to 
an ac magnetic fleld 6oe~*^* parallel to the length {x- 
direction) of the sample. Under such conditions, induced 
currents and electric flelds, J{z) and e(2;)e~*^', are along 
the y-direction, while vortices oscillate in the xz-planes. 
For low exciting flelds {bo ~ 1 /xT) , vortex displacements 
u{z)^ 1 A are very small compared with the vortex spac- 
ing a ~ n-i (- 1000 A, for B - 0.1 T) (Fig. |l|a). The 
electric fleld Cq at the surface z — d, cq — e{d) — ~e{~d), 
is measured by means of a pick-up wound coil. The main 



experimental difficulty in such measurements is to en- 
sure a precise calibration of the phase Lp of cq (within 
better than 0.5° at 100 kHz). Thus we get the sur- 
face impedance of the slab, defined as the ratio /ioCo/foo- 
Putting 

^ = ^ = y- = \'+^\" = Ke^^ , (1) 

the ac response will be conveniently expressed in terms 
of the complex penetration length A* Q. As easily seen, 
26oA (a factor 2 for two faces) represents the amplitude 
of the ac magnetic flux penetrating the slab per unit 
length along y. The length A" measures the dissipation, 
as A"/A = R/\Z\ is the sine of the loss angle tp. 

The analysis of the ideal response, though it is not ob- 
served (unless VL ^ ^d), is an important step in our argu- 
ment. A perfect slab would behave like a linear continu- 
ous medium, of resistivity pf and permeability /i = /ir/^o 
(0 < < 1). Here pf ~ pnB/Bc2 is the flux-flow resis- 
tivity, and /i is the effective "diamagnetic permeability" 
of the mixed state |^,^; Hr increases steadily with the 
vortex density and rapidly approaches unity (typically 
fir > 0.9 for B > 0.2Bc2)- In the absence of pinning, an 
electromagnetic wave, b oc e='=*'^i^e~*^', e — ^flb/ki, can 
propagate in the bulk, according to the simple equation 
of dispersion kf = ipVt/pj — 2i/5'^j Q]. The wave fleld 
b{z) would be accordingly: 610 cosh(ifci2:)/ cosh(ifcid), 
where 610 = Prbo if one makes allowance for a sur- 
face screening by diamagnetic currents. This leads to 
•^i*dcai — AirA/ tanh(c?/A/), where A/ = (1 +i)5f/2. As- 
suming /.tr ~ 1, this (undisputed) result involves all fea- 
tures of a normal skin effect. In the so-deflned thick 
limit (thick slabs and/or high frequencies), say d ^ 2(5/, 
^i*deai = MrA/ — A/, SO that A' = A" — 5f/2. In the thin 
limit, say d Sf, X^^^^i = prd ~ d, which means perfect 
transparency. 

As shown in Fig. ^, the actual response is quite dif- 
ferent: after a low-frequency plateau. A* ~ A'(0) (A" ~ 
0,ip ~ 0), the loss-angle increases with frequency, so that 
the ideal skin effect (A* ~ = 7r/4) is recovered 

beyond some depinning frequency fid] this can be pre- 
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FIG. 1. a) Vortex lines u{z) (full lines), and field profiles 
b{z) (dashed lines) near one face of a thick slab are sketched 
with arbitrary units; for clarity, the actual length scaling, 
M <C a <C Sf, is not preserved. For an ideal surface, vortices 
end normal to the plane boundary, the weight of the fc2-mode 
is lowered, and a large normal-like skin effect is observed. 
For a real rough surface, a vortex-slippage effect at the surface 
induces a relatively strong bending of vortices over the depth 
Xv ~ a, so that non-dissipative currents associated with this 
vortex curvature [^,0 are greatly screening the exciting field, 
b) An isolated superfiuid vortex in a rotating box of helium 
II terminates at a wall asperity. If it is acted on in the bulk, 
the vortex line bends near the wall so as to keep on ending 
normal to the surface, making thus an angle with the mean 
smooth surface. 



cisely defined as the mid-frequency for which (p = tt/8. 
Note, in passing, that the observation of a linear re- 
sponse is not consistent with predictions of a naive crit- 
ical state model: A critical-current density as small as 
Jc lOA/cm^ should restrict the penetration of fields 
&o ^ 1 A'T to depths A ~ bo/ ^qJc ^ 1 /im, which is much 
smaller than observed, and seeing that A cx bo, no lin- 
ear regime could exist at all. The linear skin effect over 
depths A ^ 100 /im was first reported by P. Alais and 
Y. Simon [Q| , and then misinterpreted by considering the 
possibility of thermally activated vortex motion. Soon af- 
ter, Campbell suggested that the linear signal was due to 
small reversible oscillations of vortices in their pinning 
potential wells if a pinning restoring force —nKu 
(per unit volume) is introduced in the equation of vor- 
tex motion, the propagation of the fci-mode is greatly 
altered. At low frequencies, it becomes a non-dissipative 
evanescent mode decaying exponentially in the sample 
over a small depth A(7 ~ {B(po/ fJ-aK)^ ^ 1-100 /im. Here 
na = PfK/Bifio {Sfi^d) - XcV2,ip^ tt/S) 0. Assum- 
ing /i,. ~ 1 and Ac ^ d, the Campbell expression for A* 
reads 



- J- - ± ^ J- 

1 ^ \*2 ^ \2 + \2 
A A^ A(j 



{d > 2Ac 



(2) 



With Xj = 1(5^/2, Eq. (^ accounts for the low-frequency 
plateau and the related order of magnitude of fl^- Oth- 
erwise, no satisfactory fits of both A'(il) and X"{fl) can 
be obtained from Eq. (^, as shown in Figs. || and ||. 
In spite of recent attempts to improve the treatment of 



bulk pinning, the same difficulties are encountered in fit- 
ting R{f) in YBaCuO Note in this respect that the 
inclusion of thermal fiux-creep effects Q may enhance 
the dissipation in an intermediate range of frequencies, 
as required (Fig. it should be emphasized, however, 
that flux-creep models [0,0 predicts an unobserved di- 
vergence: A' = A" oc f7-2, as 17 ^ §. 

The model of the critical state based on the Mathieu- 
Simon (MS) continuum theory of the mixed state [p|,p^, 
has prompted us to an alternative interpretation. We 
briefly recall the points of importance in the MS theory: 
i) each vortex line (unit vector u) must terminate nor- 
mal to the surface {u x n — 0); whence the leading part 
of the boundary conditions (rough or smooth surface) in 
any problem of equilibrium or motion of vortices, ii) Vor- 
tex lines are not always field lines, so that the vortex field 
uj = ntpou and the mesoscopic field B must be regarded 
as two locally independent variables. The conjugate vari- 
able of a;, £ = e{u!, T) v, appears as a local line tension 
ipo£ (J/m) in the MS equation for vortex equilibrium or 
non-dissipative motion Jg + curie = 0. iii) The classi- 
cal picture of a local diamagnetism is misleading j^. A 
diamagnetic current, just like a subcritical transport cur- 
rent, is a true non-dissipative supercurrent Js {— —curie) 
flowing near the surface over a small vortex-state pene- 
tration depth Ay (< Ao, the zero- field London depth). 
Any deviation u> — B also heals beyond the depth Ay, 
so that uj = B in the bulk sample. Although the mean 
magnetic-moment density of a perfect body turns out to 
be — e, the quantity — e has not the primary physical 
meaning of a local magnetization, nor fj,r, conveniently 
defined as the ratio uj/{uj + /ioe), that of a true local 
permeability 

Let us return to the ac response of the slab in the stan- 
dard geometry of Fig. ^ and suppose that bulk motions 
are unrestrained (no bulk pinning). As pointed out by 
Sonin et al. Q , the distinction between us and B implies 
additional degrees of freedom, and a second A;2-niodc ap- 
pears besides the classical fci-mode (fci — ±i/A/); this 
is a London-like non-dispersive evanescent mode, which 
dies off over the depth Ay (fc2 — ±i/Ay) |l^. Note that, 
in practice. Ay ^ (5/, A and d. In a pure fc2-mode vortex 
and field lines bend in opposite directions, whereas they 
coincide in a pure fci-mode. More explicitely 11|, 
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where Vx — du/dz {vix — ikiuix, V2x — ik2U2x) and 
b/B measure the slopes of vortex and field lines respec- 
tively. Then the response b{z) will be that combination 
of the modes, 61 -I- 62, which satisfies the field continuity, 
^0 = ^10 + ^20 ; well as the correct boundary condition 
for vortex lines. The latter condition will determine the 
relative weights /3i = 610/60 and (32 = 620/60 of the modes 
(/3i + P2 — 1), and, therefore, the effective penetration 
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FIG. 2. The frequency dependence of the effective penetra- 
tion depth A* = A' + in the thick limit (d > 2J/). Exper- 
imental data: (O) Pb.g2In.i8 (2d = 1.26 mm, pj = 4.8/iSl.cm, 
nd/27r = 6 MHz). (•) pure niobium (2d = 0.85mm, 
Pf = 4.3nf2.cm, Q,d/2-K — 6kHz). Full lines are theoretical 
fits using Eq.(^, where Ag is the only adjustable parameter; 
the flux-flow resistivity pf (then 5f or A/) is measured from 
the dc voltage-current characteristics. According to Eq.(P) 
the universal Argand diagram of A*(/)/A*(0) is the quarter 
circle shown in the inset: [l -I- ($1/0^)6"*'^''*] ^. For compar- 
ison the dashed line is the diagram predicted by the Campbell 
equation (0). 



depth A* = f3iXf + /32Ay ~ /3iA/ (since Ay <C(5/ and A). 
For an ideal surface, the vortex boundary condition is 
clearly i/^ — (point i) above); then using Eq. (||), 
we just recover the simple classical-diamagnetism result, 
that is Pi = fir- Now, the point is that the surface rough- 
ness may considerably change this boundary condition, 
so as to enhance the weight of the second mode (Fig. [l^) . 
Thus we argue that small effective skin depths at low 
frequencies should not result from restricting the pen- 
etration of the /ci-mode, but from its amplitude being 
reduced due to the screening effect of the second mode. 

According to the MS model §,0,||jl3|, if the sur- 
face has irregularities on the scale of a, vortex lines can 
bend over a depth Ay, making thus an angle a with 
the mean smoothed surface z — d. On the average, and 
in any direction x, a should not exceed a critical value 
1-10° : {y^),=d < iy^c= sin a, ^ IQ-^-iQ-i |,|lo|. 
As stated above, superficial non-dissipative supercurrents 
(Js = —curie) can result from such distortions of the 
vortex array; integrating —curie, the net current den- 
sity in the y-direction is Zj,(A/m) = (e^;) = e{vx)^ where 
e :^ B {I — fir) / fir ■ A dc subcritical current can be re- 
garded as a frozen pure A:2-mode in the limit ^ 0. 
To the maximum, %y = % = —VxcB{l — fir)/ fir is the 
surface critical-current density. Starting from an equilib- 
rium, where {vx) = 0, a shift of the bulk vortex array 
is expected to induce a vortex curvature in the opposite 
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FIG. 3. The complex penetration depth A* = Ae*"^ of 
the Pbin slab referred to in Fig. |2| has been measured as a 
function of the magnetic fleld B, for three values of the driving 
frequency. Experimental data are plotted as the sine of the 
loss angle ip as function of the ratio A/prSf, so as to verify the 
relation A — prSf{^l, B) simp (straight line) resulting from 
Eq.(||). The limit sin 93 = 1/^2 corresponds to the normal 
state or to the depinned vortex array. The Campbell equation 
(^) leads to much smaller loss angles: A^ = 1(5^ sin 2ip (dashed 
line). 



direction, = /(ubuik) with v^c = /(u~a). Perhaps 
such a vortex slippage {vx < if u > 0) is more intuitive 
when dealing with superfluid vortices in a rotating box 
(see Fig. |l|b). In helium, bulk pinning does not exist, and 
only asperities of the walls can pin vortices |14|. We are 
just extending this idea to collective motions of a vortex 
array along a rough surface of superconductors. 

Linearizing Vx = /(ubuik) for small displacements we 
can write Vx — — Wbuik/^, where I is a real length char- 
acterizing the surface roughness: I = 00 would corre- 
spond to an ideal surface; in practice we expect that 
I ^ a/i^c ^ 0.1 — 10 fim. As far as Ay ^ 5f, this condi- 
tion applies quasistatically in the ac response by taking 
Mbuik = Mio, so that the vortex-slippage condition reads 



^10 

I 



(4) 



Now, substituting Eq. (Q) for i/^, = in the above calcu- 
lation of Pi and P2, and considering the thick limit, we 
obtain: 
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H , {thick limit, d}^ 26 f) (5) 
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where A5 = //^^/(l — /i^) ^ Ba/fiQ^c is the real limit of 
A* ~ PiXf as ^ (A/ = (1 -h i)5f/2). Note that 
setting /i^ = 1 from the beginning would lead wrongly 
to A* = A/, that is the ideal response. While giv- 
ing the same low and high-frequency limits as Eq. (m), 
A*(0) = As (real) and A* (00) = /irA/ (depinning), Eq. ^) 



3 




FIG. 4. The microwave surface resistance R of an YBaCuO 
film (sample No. 2 of Ref. Q) vs frequency plotted in a log-log 
scale. (•) are experimental data taken from Fig. 6 of Ref. 
The full line is a fit using our Eq. and taking As and p/ 
(or Xf) as two adjustable parameters (we find As = 0.07/^m, 
pf = OAfiflcm, fld/2-K = 100 GHz). This fit is very close to 
the empirical power-law R ~ f^'^'' proposed by the authors. 
The dashed line shows the best fit obtained in Ref. from 
the Coffey-Clem fiux creep-model |^ 

remarkably fits experimental data in the intermediate 
range fid (Figs. ^ and ||). According to Eq. (||), the 
graph A*(ri) in the Argand diagram must be a quarter 
circle (Fig. ||). This universal behaviour should be eas- 
ily tested in any case, irrespective of any adjustable or 
available parameters, only providing that the thick Umit 
is achieved. 

From data taken below the thick limit (not shown 
in the figures), we retain an important result. When 
2d is decreased under conditions where Ac or A5 ^ (5/ 
(fi ^ Od), we observe that the ac response becomes 
thickness-dependent as soon as d < 2Sf, as predicted: 
just substitute A/tanh(d/A/) for A/ in Eq. (^). In par- 
ticular, the loss angle is significantly smaller than stated 
by Eq. (^) . According to the one-mode Campbell model, 
size effects should arise for much thinner slabs such as 
d < Ac- Note in this respect that our "pinning length" 
As, contrary to Ac, does not represent an actual penetra- 
tion depth. The mere observation that size effects arise 
for d ^ 5f, not d ~ Ac, attests that the bulk /ci-mode 
propagates freely, and reveals the need for a two-mode 
electrodynamics. 

In conclusion, the MS model of the critical state 
completed by the vortex-slippage condition (^) accounts 
quantitatively for the surface impedance of a variety of 
conventional samples, which all have standard critical- 
current densities (Zc ~ 10 A/cm): polycristalline lead- 
indium slabs and single-crystal slabs of pure niobium 
(Figs. |2| and For the application to the case of 

YBaCuO at / ~ 10 GHz, our derivation of Eq. (||) has to 



be reexamined, especially because of the anisotropy and 
high-frequency correcting terms in the dispersion equa- 
tion for the two-modes ||ll|. Nevertheless, it is worth 
noting that Eq. (||) may account for the observed /- 
dependence of R in YBaCuO from 1 to 20 GHz (Fig. |) 
Q. These results support the argument, developed in 
previous works, ^,|l^,|l2|,|l3| that bulk pinning is abso- 
lutely ineffective in a large class of "soft" materials (de- 
void of strong bulk inhomogeneities). Contrary to the 
common idea that any crystal defect may be a pinning 
center, we are led to the conclusion that a normally ho- 
mogeneous sample in the mixed state rather imitates the 
behaviour of a superfiuid vortex array enclosed in a rough 
box. 
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